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Abstract:  Laser  enables  the  achievement  of  superb  interfacial  characteristics  between 
electrode  and  semiconducting  material  contact  surface  and  is  also  useful  for  a  reduction  in 
contact  resistance.  The  irradiation  of  a  pulsed  laser  with  high  energy  density  and  short 
wavelength  onto  the  electrodes  leads  to  the  thermal  annealing  effect  at  the  locally  confined 
small  area  that  needs  high  temperature  without  extreme  thermal  damage.  Lfnlike  pulsed  laser 
annealing,  conventional  thermal  annealing  affects  the  entire  panel  including  unwanted  areas 
where  the  annealing  process  should  be  excluded. 

2D  transition  metal  dichalcogenides  (TMDCs)  with  the  formula  of  MX2  (M  =  Mo,  W; 
X  =  S,  Se,  Te)  composed  of  vertically  stacked  layers  held  together  by  van  der  Waals 
interactions,  have  shown  potential  to  complement  or  overcome  drawbacks  of  conventional 
opto-  and  electronic  device  applications  (e.g.  transistors,  photodetectors,  electroluminescent 
devices  and  sensors),  because  of  their  unique  optical  and  electrical  properties:  relatively  large 
bandgap  (1.2- 1.9  eV),  high  mobility  at  room  temperature  (up  to  -500  cm^V'^s  '),  low 
subthreshold  swing  (SS,  -70  mV  decade'^),  high  on/off  ratios  (-10^),  and  potential 
amenability  to  low-cost  fabrication  sequence.  Tuning  their  properties  by  modulating  the  free 
carrier  type  and  density  and  composition  can  offer  an  exciting  new  pathway  to  various 
practical  nanoscale  electronics.  However,  there  are  several  key  challenges  that  lie  ahead 
before  the  transistors  based  on  2D  layered  semiconductors  can  be  integrated  into 
flexible/stretchable  circuits.  The  main  issues  include  instability  of  flexible  structures  with 
respect  to  various  mechanical  stresses,  inaccessibility  of  large  area  growth/coverage  of 
layered  semiconductors,  and  tunability  of  TMDC-based  layered  semiconductors  for 
ultra-sensitive  sensors  of  high  specificity.  These  challenges  give  rise  to  the  focus  areas 
defined  in  our  research  plan  for  this  project. 

Experimental  and  computational  studies  in  multidisciplinary  fields  of  electrical, 
mechanical  engineering,  and  materials  science  were  conducted  to  achieve  the  following 
major  goals  of  this  project: 

•  Development  of  laser  annealing  for  improving  contact  quality  to  achieve 
high-performance  devices  on  a  flexible  platform 

•  Modification  of  optical  and  electrical  properties  of  2D  materials  with  laser 
assisted  and  site  selective  doping  and  alloying 

•  Development  of  plasmonic  laser  welding  process  and  solution  based  embedded 
polyimide  (PI)  structure  for  making  highly  robust  flexible/wearable  devices 

•  Stress-strain  analysis  for  optimizing  the  structure  design  (e.g.  radius  of  curvature, 
thickness  of  the  structure,  etc.) 

•  Large-area  and  highly  crystalline  CVD-grown  2D  films 

•  Optimization  of  device  designs  to  achieve  flexible  biosensors  with  2D  materials 
for  highly  sensitive  and  selective  biological  detection 


Through  the  project,  aforementioned  results  demonstrate  that  pulsed  laser  process  is 
an  attractive  technology  for  realizing  high  performance  flexible  M0S2  TFTs  in  analog/digital 
integrated  circuits;  they  are  published  in  peer-reviewed  journals  (more  than  17  papers)  and 
are  presented  at  conferences  (more  than  10  times).  In  addition,  we  are  continuing  our  research 
on  flexible/wearable  electronics  and  their  applications,  especially  for  sensors. 
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Introduction:  Flexible  and  stretchable  electronics  represent  a  critical  frontier  in  the 
transformation  of  rigid,  table-top  micro/nano  electronics  into  portable,  wearable  systems  that 
can  be  integrated  into  a  variety  of  emerging  technologies  from  sensing  and  monitoring  to 
human-inspired  applications.  Many  conventional  structures,  materials,  and  processes  are  not 
compatible  with  flexible/stretchable  device  layouts.  These  requirements  pose  significant 
challenges  that  require  a  new,  adaptive  paradigm  for  the  low-thermal  budget  (<  100  °C)  and 
functional  components  on  a  lightweight  and  inexpensive  flexible  platform.  Thin-film 
transistors  (TFTs)  are  representative  applications  in  the  field  of  flexible/stretchable 
electronics.  Currently,  users  require  new  and  more  advanced  displays  for  mobile  phones, 
portable  devices,  and  even  televisions,  with  features  including  ultra-high  resolution,  higher 
frame  rates  (high  driving  speed),  low  power  consumption,  and  larger  size.  For  realizing  high 
performance  of  TFTs,  we  should  consider  metal-semiconductor  junctions  that  are 
indispensable  to  the  structure  of  TFTs,  because  non-ideal  electric  contacts  can  form  Schottky 
barriers  at  the  junctions,  and  the  undesirable  contacts  can  hamper  the  inherent  electrical 
characteristics  of  M0S2.  Therefore,  to  reduce  the  contact  resistance  (Rc)  and  Schottky  barrier, 
various  efforts  (doping,  the  use  of  scandium  metal  with  a  low  work  function,  and  thermal 
annealing)  have  been  reported.  However,  there  are  challenges  and  limitations  to  moving 
forward  in  flexible  electronics  technologies.  The  chemical  doping  effect  is  gradually  reduced 
over  time.  Moreover,  scandium  metals  have  been  classified  as  a  rare  earth  element  and 
difficulties  have  been  encountered  in  preparing  them.  The  typically  used  flexible  substrates 
(e.g.  poly(ethylene  naphthalate)  (PEN),  poly(ethylene  terephthalate)  (PET),  and  polyimide 
(PI)),  have  a  low  thermal  budget  (<  200  °C).  Therefore,  high  thermal  annealing  processes 
could  lead  to  deformation  in  plastic  substrates. 

Unlike  conventional  thermal  annealing,  which  affects  the  entire  panel  including 
unwanted  areas  where  the  annealing  process  should  be  excluded,  the  irradiation  of  a  pulsed 
laser  with  high  energy  density  and  short  wavelength  onto  a  metal  electrode  leads  to  the 
thermal  annealing  effect  in  a  small,  locally  confined  area  that  requires  high  temperature 
without  extreme  thermal  damage.  Also,  a  laser  enables  the  achievement  of  superb  interfacial 
characteristics  between  metal  and  semiconducting  material  contact  surfaces,  resulting  in 
reduced  Rc  and  improved  interfacial  morphology.  Numerical  thermal  simulation  for  the 
temperature  distribution,  transmission  electron  microscopy  (TEM)  analysis,  current-voltage 
measurements,  and  contact-free  mobility  extracted  from  the  Y-function  method  (YFM) 
enable  understanding  of  the  compatibility  and  the  effects  of  pulsed  laser  annealing  process; 
the  enhanced  performance  originated  not  only  from  a  decrease  in  the  Schottky  barrier  effect 
at  the  contact,  but  also  an  improvement  of  the  channel  interface. 

In  addition,  doping  is  simple  and  versatile  methods  to  accomplish  our  goal.  However, 
previous  efforts  have  mostly  focused  on  doping  TMDCs  by  means  of  charge  transfer  from 
adsorbed  molecules,  electrostatic  or  physisorption  gating,  defect  engineering,  and 
substitutional  doping  during  the  growth.  These  methods  result  in  uniform  doping  or 
composition  without  spatial  resolution.  However,  for  nanoscale  devices  using  ultrathin 
TMDCs  as  an  ultra-thin  body,  to  minimize  random  dopant  fluctuation,  ensure  device 
performance  reproducibility,  and  scale  down  the  device  size  and  scale  up  device 
functionalities,  site-specific  doping  with  precise  carrier/composition  level  control  becomes 
essential.  Here  we  provide  a  systematic  approach  to  this  problem,  demonstrating  superior, 
in-situ  control  of  the  spatial  distribution  and  doping  level  of  TMDCs.  Our  laser  induced 
doping  process  is  spatially  local  and  selective,  and  the  doping  level  is  widely  tunable  by 
varying  the  duration  and  intensity  of  laser  irradiation.  Moreover,  laser  doping  can  be  digitally 
controlled  and  is  an  air-stable  process  of  TMDCs  for  high-quality,  high-fidelity  nanodevices. 

Many  research  groups  have  attempted  to  develop  alternative  conductive  flexible 
electrodes,  since  the  usual  metal  films  are  rather  ductile  and  require  expensive  vacuum 
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deposition  techniques.  Candidate  electrode  systems  that  have  been  studied  and  introduced  up 
to  date  include  highly  conductive  poly(3,4-  ethylenedioxythiophene):poly(4-styrenesulfonate), 
metal  grids,  graphene,  carbon  nanotubes,  nanotube-  polymer  composites,  and  silver  nanowire 
(Ag  NW)  meshes.  Among  the  suggested  materials,  2D  random  networks  of  Ag  NWs  are 
promising  because  of  their  compatibility  with  cost-effective  solution  processes  and  the 
reasonably  high  electrical  conductivity.  However,  random  network  of  Ag  NWs  may  be 
inappropriate  for  practical  applications.  First,  the  electrical  properties  must  be  preserved  upon 
mechanical  bending  or  stretching.  Next,  the  surface  roughness  of  the  entire  electrode  layer 
must  be  sufficiently  small  in  order  to  guarantee  fast  carrier  transport  along  the  overlaying 
active  layers.  Finally,  the  NWs  must  adhere  well  to  the  substrate  and  form  contacts  at  their 
intersections.  The  electrical  contact  at  the  Ag  NW  junctions  can  be  improved  by  laser 
welding,  which  is  based  on  highly  localized  heating  induced  by  plasmonic  effects. 


Experiment: 

a.  Laser  annealing 

Fig.  1(a)  shows  the  schematic  architecture  of  our  flexible  multilayer  M0S2  FETs  on  the  PEN 
substrate  (DuPont  Teijin  Eilms,  USA). 
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Figure  1:  Flexible  multilayer  M0S2  FETs  and  pulsed  laser  annealing  system,  (a)  Cross-sectional  schematic 
structure  of  multilayer  M0S2  FETs  with  a  300-nm- thick  Si02  gate  dielectric.  Picosecond  laser  irradiation  with  a 
Gaussian  beam  profile  is  focused  and  scanned  along  the  edge  of  the  source  and  drain  electrodes,  (b)  Top-view  of 
the  optical  microscope  image  for  multilayer  M0S2  FETs  after  fabricating  completely,  (c)  Schematic  view  of  the 
optical  setup  for  picosecond  laser  annealing. 
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For  a  bottom  gate  electrode,  a  100-nm-thick  indium  thin  oxide  (ITO)  was  deposited 
on  the  PEN  substrate  at  room  temperature.  Then,  dual  films  (atomic-layer-deposited  (ALD) 
AI2O3  of  50  nm  /  sputter-deposited  SiOa  of  250  nm)  were  formed  as  a  gate  dielectric  layer.  A 
M0S2  flake  with  a  thickness  of  ~  50  nm  was  mechanically  transferred  from  bulk  M0S2 
crystals  (SPI  Supplies,  USA)  on  the  deposited  dielectric  layer.  After  that,  Ti  (10  nm)  and  An 
(300  nm)  were  continuously  deposited  by  electron-beam  evaporation  at  room  temperature. 
Finally,  conventional  UV  photolithography  and  lift-off  techniques  allowed  definition  of  the 
source  and  drain  electrodes.  Fig.  1(b)  shows  a  proof-of-concept  realization  of  a  completely 
fabricated  flexible  multilayer  M0S2  FET  rolled  up  in  a  curved  shape  and  an  optical 
microscope  (BX60,  Olympus)  image  of  the  top  of  the  device  (Fig.  lb  inset).  A  yttrium 
vanadate  (Nd:YV04)  picosecond  pulsed  laser  (Newport  Corp.,  USA)  having  a  wavelength  of 
355  nm,  a  pulse  width  based  on  full  width  at  half  maximum  of  12  ps,  and  a  pulse  repetition 
rate  of  80  MHz  was  applied  for  selectively  focused  annealing  at  the  contact  regions.  Note  that 
the  shape  of  the  pulsed  laser  beam  was  Gaussian  with  a  1.5-pm  beam  diameter  at  a  1/e^  peak 
irradiance,  and  was  obtained  through  a  39x  objective  lens  with  a  numerical  aperture  (NA)  of 
0.5.  Moreover,  a  computer  simultaneously  operated  laser  power,  scan  speed,  position  of  the 
high-resolution  x-y  positioning  stage  (Aerotech,  USA),  and  optical  shutter  during  the 
annealing  process,  as  shown  in  Fig.  1(c). 

Additionally,  the  Korean  partner  fabricated  TFTs  based  on  multilayer  M0S2  using 
mechanically  exfoliated  multilayer  M0S2  flakes  on  a  highly-doped  p-type  Si  wafer  with  a 
thermally  grown  Si02  layer  (300  nm).  Electrical  contacts  (100  pm  x  100  pm)  were  patterned 
on  top  of  M0S2  flakes  using  conventional  photolithography  and  lift-off  techniques.  Ti  (100 
nm)  was  deposited  by  electron-beam  evaporation  at  room  temperature.  The  device  was  then 
annealed  at  400°C  in  a  vacuum  tube  furnace  for  2  hours  (100  seem  Ar  and  10  seem  H2).  The 
thickness  of  M0S2  channels  measured  by  atomic  force  microscope  (AFM;  Park  Systems, 
XE-100)  was  in  the  range  of  20-80  nm.  Electrical  characterizations  were  carried  out  with  I-V 
measurements  (Keithley,  Semiconductor  Characterization  System  4200-SCS)  in  vacuum. 

b.  Faser  doping 

The  schematic  diagram  of  the  laser-assisted  doping  method  shows  in  Fig.  2.  The  laser  serves 
two  major  functions  in  this  doping  method,  one  is  to  create  surface  vacancies  in  the  TMDC 
materials  and  the  other  is  to  dissociate  the  dopant  gas,  releasing  the  dopant  molecules  to 
incorporate  into  the  vacancy  sites.  The  phosphine  (PH3)  was  introduced  as  p-type  dopant 
precursor  for  ultrathin-layered  materials,  in  contrast  to  its  use  for  n-type  doping  silicon.  The 
PH3/Si  system  has  been  extensively  studied  both  experimentally  and  through  theoretical 
modeling.  On  the  basis  of  this  insight,  and  despite  the  structural  difference  between  silicon 
and  TMDCs,  it  is  possible  to  adopt  the  vacancy  mechanism  in  this  work.  In  terms  of  vacancy 
formation,  it  is  noted  that  laser  power  just  below  dissociation  of  the  TMDC  materials,  is 
enough  for  breaking  the  PH3  (its  dissociation  temperature  is  685  K)  since  the  TMDC  material 
dissociation  temperature  is  in  the  range  of  1200-1400  K. 


Figure  2:  Schematic  diagram  of  the  laser  assisted  doping  method.  A  mechanically  exfoliated  M0S2  flake  is 
transferred  onto  the  300nm  Si02  on  Si  substrate  and  Phosphine  (PH3)  as  a  dopant  gas 
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c.  Laser  welding 

To  make  a  Ag  NWs  network,  a  vacuum  transfer  method  was  applied  because  the  vacuum 
suction  pressure  can  distribute  the  normal  force  uniformly  on  the  substrate.  Firstly,  a  Teflon 
filter  (SterliTech,  0.2  pm  pore  size  and  47  mm  in  diameter)  was  placed  at  the  neck  of  a  side 
arm  flask  connected  to  a  vacuum  pump  and  the  glass  cap  was  put  on  the  Teflon  filter  for 
filtration.  Then,  25  pL  of  Ag  NW  solution  (Seashell  Technology,  with  25  +  10  pm  length  and 
30+10  nm  diameter)  dispersed  in  sufficient  isopropyl  alcohol  (as  a  carrier  fluid)  was  poured 
into  the  cap.  Note  that  Ag  concentration  was  around  5.3  mg  ruL”'.  Only  Ag  NWs  were 
filtered  on  the  Teflon  filter.  After  drying  for  3  hours,  the  filtered  Ag  NWs  were  directly 
transferred  from  the  filter  to  the  target  substrate  (Glass- 1  in  Fig.  3)  with  a  uniform  normal 
pressure  for  1-2  hours.  Owing  to  very  low  surface  energy  of  Teflon,  most  of  Ag  NWs  were 
transferred  easily  and  a  uniform  2D  Ag  NWs  percolation  network  was  formed  on  the  target 
substrate.  The  as-prepared  Ag  NWs  had  physical  contacts,  which  is  susceptible  to  the  shear 
force  and  can  result  in  high  contact  resistance.  Therefore,  to  enhance  the  electrical  property 
by  forming  good  contacts,  a  laser  beam  with  532  nm  wavelength  (Solo  I  Nd:YAG,  New 
Wave)  was  irradiated  onto  the  Ag  NWs  network.  Also,  through  an  objective  lens,  the 
effective  area  was  scaled  up  to  -cm  scale  by  an  out-of-focal-plane  process  with  a  maximum 
laser  energy  of  20  mJ. 

For  a  flexible  substrate,  PI  solution  was  spin-coated  on  the  surface  of  the  laser-welded 
Ag  NWs  at  3,000  rpm  for  30  s  and  thermally  solidified  at  60  °C  for  30  min.  The  Ag 
NWs-embedded  PI  layer  was  mechanically  detached  from  the  Glass- 1,  and  attached  on  the 
carrier  substrate  (Glass-2)  upside  down  with  3M  tapes  at  each  corner.  The  organic  GI  was 
formed  by  spin-coating  SU-8  2000.5  (Microchem®)  onto  the  Ag  NWs  gate  electrode  (3,000 
rpm  for  30  s),  followed  by  baking  at  110  °C  for  5  min  and  130  °C  for  5  min.  Then,  a 
30-nm-thick  amorphous  AI2O3  was  deposited  by  ALD  (Lucida  DlOO,  South  Korea)  with 
trimethylaluminum  as  an  A1  precursor  and  H2O  as  an  oxidant.  Mechanically  exfoliated 
multilayered  M0S2  flakes  were  transferred  onto  the  AI2O3  film.  After  E-beam  evaporation  of 
a  Ti/Au  (20/100  nm)  layer,  source  (drain)  electrodes  were  fabricated  by  conventional 
photolithography  and  etching  techniques.  The  as-fabricated  multilayer  M0S2  TFTs  on 
flexible  PI  substrate  was  post-annealed  at  150  °C  for  2  hours  in  vacuum  oven  (<  10'^  torr)  in 
order  to  reduce  contact  resistance  between  the  active  layer  and  the  S/D  electrodes. 
Transparent  PET  film  with  a  thickness  of  100  pm,  as  a  protective  as  well  as  supportive  layer 
for  the  flexible  M0S2  TETs,  was  laminated  beneath  the  PI  layer  with  the  aid  of  an  acrylate 
resin-based  pressure  sensitive  adhesive  after  the  separation  of  the  flexible  M0S2  TETs  from 
the  carrier  substrate. 


Ag  NWs 


Vacuum  Transfer  of  Ag  NWs~~j 


Separating  TFT  from  Glass-2 1 
Lamination  of  Protective  Film 


Figure  3:  Flexible  multilayer  M0S2  TFTs.  Schematic  flow  chart  of  the  fabrication  procedure  consisting  of  the 
laser  welded/embedded  2D  Ag  NWs  network,  solution-based  thin  PI,  and  organic/inorganic  hybrid  GI  for  a 
highly  flexible  TFT  platform. 
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Results  and  Discussion: 

a.  Laser  annealing 

i.  Thermal  analysis  for  pulsed  laser  annealing  system 
Profound  understanding  of  ultrashort  laser  pulse  induced  temperature  field  requires  detailed 
knowledge  of  the  fundamental  interactions  between  laser  light  (photons)  and  matter  (electron 
and  phonon  systems).  In  metals,  initially  electrons  primarily  absorb  the  laser  photon  energy. 
In  the  Au  film  with  the  thickness  of  400  nm,  the  optical  absorption  depth  is  about  16.3  nm. 
The  electron  system  temperature  rises  instantaneously  in  response  to  the  applied  photon  flux. 
Thermalized  hot  electrons  diffuse  into  the  material  and  transfer  their  energy  to  the  lattice 
(within  lO  '^  ~  10  '°  s)  until  the  electron  and  lattice  temperatures  equilibrate  through 
electron-phonon  coupling.  This  equilibration  occurs  within  the  pulse  for  pulse  lengths  in  the 
range  10  to  60  ps.  In  the  gold  system,  this  thermal  equilibrium  time  of  the  electron  and  lattice 
is  about  40  ps.  Therefore,  the  temperatures  of  electron  and  lattice  are  different  during  the 
single  pulse  duration  (12  ps)  applied  in  our  experiment.  Consequently,  the  hot  electron 
diffusion  length  is  deeper  than  both  the  optical  absorption  depth  and  the  heat  diffusion 
length  8~d4adt  (where  d  is  the  depth  of  heat  penetration,  a  is  the  diffusivity  of  materials, 
and  dt  is  the  pulse  duration).  In  order  to  obtain  a  more  accurate  predictive  solution,  the 
two-temperature  model  (TTM)  could  be  applied  in  this  transient  non-equilibrium  case. 
However,  the  TTM  model  that  is  by  itself  an  approximation  of  the  Boltzmann  transport 
equation  has  been  widely  used  to  understand  the  origin  of  ultrashort  electron  dynamics  under 
single  laser  pulse  excitation  of  width  lO  '^  ~  lO  '^  s.  Prediction  of  the  transient  temperature 
distribution  induced  by  high  repetition  rate  pulsed  laser  over  a  long  time  scale  as  in  our 
system  by  using  the  TTM  poses  a  significant  challenge  due  to  the  required  prohibitively  long 
computing  time.  In  our  80  MHz  repetition  rate  case,  we  have  to  compute  the  TTM  system  of 
equations  over  80,000,000  pulses  in  order  to  analyze  only  one  second.  Although  the 
amplitude  of  the  ps  duration  temperature  spikes  would  be  smaller  by  the  TTM  model,  the 
heat  accumulation  that  is  driven  by  thermal  diffusion  would  not  differ  from  the  presented 
thermal  model  prediction. 

For  a  picosecond  pulsed  laser  source  applied  under  our  experimental 
conditions,  we  investigated  temperature  distribution  during  laser  annealing  process  by 
3-dimensional  (3D)  finite-difference  methods  (FDM)  using  COMSOL  Multiphysics. 
The  equilibrium  temperature  approximation  is  adopted  since  it  reduces  the  required 
computational  cost.  As  previously  noted,  the  assumption  of  equilibrium  temperature 
would  overestimate  the  transient  temperature.  However,  the  approximate  thermal 
solution  shows  good  agreement  with  the  results  of  real  laser  damage  test.  To  calculate 
the  approximated  temperature  distribution  inside  laser-heated  Au  electrode,  the 
following  heat  flow  equation  is  used: 

pC/^  =  V{k(T)VT)  +  Qix.y.zy  (1) 


where  p  is  the  density,  Q,  is  the  specific  heat,  k  is  the  thermal  conductivity,  and  Q  is 
the  heat  source  term.  In  the  volumetric  heating  with  a  Gaussian  laser  beam  case,  the 
heat  source  term  is  expressed  as: 


Qix,y,z)  = 


Qo(i-Rc)- 


2jtoa^, 


(2) 


where  Qo,  Rc,  Ac,  and  Oxj  indicate  the  total  input  power,  the  reflection  coefficient,  the 
absorption  coefficient,  and  pulse  x,y  standard  deviation,  respectively.  Note  that,  for 
the  accuracy  of  the  analysis,  temperature  dependent  thermophysical  properties  of  Au 
electrode  used  as  absorption  layer  were  applied.  However,  other  material  properties. 
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Rc  and  Ac,  are  assumed  to  be  constants.  In  addition,  the  thermal  properties  of 
multilayer  M0S2  were  included  in  this  analysis  by  using  the  experimental  values. 

The  top  surface  of  An  electrode  is  aligned  with  z  =  0  and  the  Gaussian  laser 
beam  moves  over  the  surface  at  a  specific  speed,  10  pm  s  '  along  a  prescribed  path 
(Fig.  4(a))  because  charge  injection  mainly  occurs  from  the  edge  of  the  metal 
electrodes  into  the  semiconducting  channels.  In  addition,  the  temperature  generated 
by  the  laser  annealing  did  not  cool  down  to  below  ambient  temperature  before  the 
next  pulse  arrives  due  to  the  high  repetition  rate  (80  MHz).  This  heat  accumulation 
effect  can  induce  changes  in  the  characteristics  of  M0S2  FETs  and  have  made  a  pulsed 
laser  usable  for  annealing  flexible  electronics.  Fig.  4(b)  shows  the  optical  images  of 
laser-irradiated  An  contacts  on  plastic  substrate.  The  laser  power  over  40  mW 
severely  damaged  the  An  electrode  because  generated  heat  could  exceed  the  melting 
temperature  of  An  (1330  K).  To  predict  the  generated  heat  and  temperature 
distribution  inside  the  sample,  we  computed  the  temperature  with  respect  to  the 
scanned  time  at  particular  points.  Note  that  we  calculate  the  temperature  distribution 
for  average  laser  powers  of  18.3  mW  (orange  region)  and  43.0  mW  (blue  region); 
18.3mW  makes  the  best  reduction  of  contact  resistance  and  43.0  mW  is  the  initial 
power  of  damage  based  on  real  experimental  results  (Fig.  4(b)).  While  the  moving 
pulsed  laser  beam  with  the  speed  of  10  pm  s  '  directly  irradiated  on  the  An  surface, 
the  effect  of  the  high  frequency  transients  by  an  ultrashort  pulsed  laser  lasted  till  the 
time  of  0.2  s;  the  temperature  distribution  then  consists  of  high-frequency  and 
picosecond-duration  spikes  superposed  on  a  continuous  base  as  shown  in  the  inset  of 
Fig.  4(c).  In  agreement  with  the  real  test  of  laser  irradiation  (Fig.  4(b)),  for  the  laser 
power  of  43.0  mW,  the  temperature  exceeded  the  melting  point  of  An  (1330  K)  as 
indicated  by  the  orange  region  in  Fig.  4(c).  Such  a  high  temperature  may  inflict 
thermal  damage  on  the  An  contacts.  However,  the  induced  temperature  by  the  laser 
power  of  18.3  mW  is  533  K  (blue  region  in  Fig.  4(c)).  This  is  a  comparable 
temperature  with  those  of  thermal  annealing  used  in  previous  studies  to  improve 
device  performance,  as  well  as  with  our  experimental  results  detailed  in  subsequent 
sections.  Also,  generated  heat  by  laser  absorption  at  the  top  surface  could  affect  the 
interfaces  by  heat  conduction  as  well  as  the  contact.  Therefore,  we  expand  thermal 
analysis  into  the  interfaces.  Note  that  interested  interfaces  are  between  channel  and 
dielectric  layer  (point  2)  and  between  dielectric  and  substrate  layer  (point  3)  in  Fig 
4(d).  In  Fig.  4(c),  we  already  analyzed  point  1.  As  shown  in  Fig  4(d),  at  point  2,  the 
reached  heat  around  480  K  (-210  °C)  is  not  low  to  ignore  because  the  induced 
temperature  is  high  enough  for  thermal  annealing,  according  to  other  reports. 

Unfortunately,  we  cannot  obtain  the  result  measured  directly  to  confirm  the 
effect  on  the  interfaces.  However,  we’ve  tried  to  verify  the  impact  on  the  interfaces 
indirectly  through  the  new  approach,  YFM  (see  the  next  session).  Point  3  shows  that 
the  reached  heat,  365  K  (-90  °C)  can  be  below  the  thermal  resistance  of  PEN 
substrate  (-170  °C).  Therefore,  we  can  know  and  expect  that  there  is  no  thermal 
damage  after  laser  process.  We  already  proved  it  through  the  laser  experiment  (Eig 
4(b)). 
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Figure  4:  Numerical  thermal  analysis  of  temperature  distribution  under  a  moving 
picosecond-pulsed  laser  source,  (a)  Spatial  schematic  images  of  laser-irradiated  Au  contacts  at  a 
speed  of  10  pm/s  on  a  flexible  PEN  substrate  and  three  points  of  interest:  point  1  (the  interface  between 
the  contact  and  M0S2),  point  2  (the  interface  between  M0S2  and  the  dielectric  layer),  and  point  3  (the 
interface  between  the  dielectric  layer  and  the  substrate),  (b)  the  optical  images  of  the  front  and  back 
side  of  laser-irradiated  Au  contacts  on  a  flexible  plastic  substrate.  Au  contacts  irradiated  with  laser 
power  >  40  mW  show  severe  damage,  (c,  d)  predicted  temperature  distributions  versus  exposure  time 
at  the  three  specific  points  when  a  picosecond  pulsed  laser  was  applied  and  scanned  onto  the  top 
surface  (Au)  at  ambient  conditions. 


a.  TEM  analysis:  Schottky  barrier  at  contacts 

Fig.  5  shows  the  cross-sectional  TEM  image.  Through  the  TEM  images  before  and 
after  laser  annealing,  there  are  several  possible  explanations  for  the  effects  of  laser 
process  on  the  flexible  multilayer  M0S2  EET.  The  first  possibility  is  removing 
unintentional  interfacial  layers  that  introduce  interface  traps.  The  second  possible 
explanation  is  a  mixture  layer.  It  is  not  clear  which  mechanism  makes  a  mixture  layer, 
however,  the  mixture  layer  is  conspicuously  formed  after  laser  annealing  process  as 
shown  in  Eig.  5.  This  observation  is  consistent  with  the  reports  that  Ti  forms  a 
mixture  layer  with  M0S2  by  placing  Ti  atoms  between  the  neighboring  S  planes  of 
M0S2  and  that  thermal  annealing  at  300  °C  induces  solid-state  diffusion  in  Ti  thin 
film. 
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Figure  5:  Effect  of  laser  annealing  on  metal-MoS2  junction  before  and  after  laser  annealing. 

Cross-sectional  TEM  image  showing  the  mixture  layer  at  the  interface  after  laser  annealing. 


Hi.  Electrical  characteristics 

Fig.  6(a)  shows  the  log-  and  linear-scale  drain  current  vs.  gate-source  voltage  (Ids-Vgs) 
characteristics  for  the  same  M0S2  TFTs,  measured  before  (black  circles)  and  after 
laser  annealing  (red  circles)  with  a  power  of  18.3  mW).  Following  selective  laser 
annealing  on  contacts  in  M0S2  TFTs,  transistor  performance  metrics,  including 
on/off-current  ratio  (lon/Ioff),  peak  field-effect  mobility  (pfe  =  Lgm/(WCoxVds)), 
subthreshold  swing  (SS),  output  resistance  (Ro),  and  voltage-gain  for  transistor  (Ay), 
were  improved.  The  Wloff  after  laser  annealing  was  observed  as  >  10^  without  a 
change  in  loff.  pfe  extracted  from  Fig.  6(a)  increased  almost  two-fold  after  laser 
annealing  (from  24.84  cm^  V'^s'^  to  44.84  cm^  V's’').  Here,  L,  gm,  W,  Cox  are  channel 
length  (L  =  7  pm),  transconductance,  channel  width  (W  =  7.75  pm),  and  oxide 
capacitance,  respectively.  The  value  of  pfe  observed  in  our  M0S2  TFTs  exceeds  those 
reported  in  previous  M0S2  TFTs  on  flexible  substrates  (4.7-24  cm^  V'^s'^)  as  well  as 
in  organic  TFTs  or  a-Si  TFTs  on  flexible  plastic  substrates  (<  1  cm^  V's  '). 
Statistically,  pfe  of  the  tested  M0S2  TFTs  increased  from  -15+10  cm^  V's'^  to  -35+10 
cm^  V'^s  '  after  laser  annealing.  No  dependence  on  channel  thickness  was  observed 
within  the  given  experimental  conditions.  Fig.  6(b)  shows  the  drain  current  as  a 
function  of  the  drain  voltage  (Ids-Vds)  with  magnification  for  the  low  Vds  range  (from 
0  V  to  0.4  V)  before  (open  circles)  and  after  laser  annealing  (solid  circles).  At  the 
same  Vgs-Vth,  the  laser  annealing  results  in  ohmic  behavior  in  linear  regime  and 
higher  current  saturation.  Following  laser  annealing,  Ro  significantly  increased  to  0.42 
MQ  at  Vgs-Vth=  20  V,  which  is  three  times  higher  than  values  for  as-fabricated 
devices.  The  most  prominent  effect  of  laser  annealing  is  the  self-gain  (Ay  =  gmRo)  in 
M0S2  TFTs,  which  exhibits  six-fold  increase  due  to  enhancement  of  gm  and  Rq. 
Furthermore,  as  shown  in  the  Fig.  6(c),  we  can  observe  the  non  linear  output  curve 
due  to  the  Schottky  barrier  before  laser  process.  However,  after  laser  annealing,  the 
curve  exhibited  linear  characteristics.  It  could  be  possible  to  reduce  the  barrier  effect. 
For  supporting  the  results,  we  already  examined  the  effects  of  the  laser  annealing 
through  the  cross-sectional  TEM  image;  the  reduction  of  tunneling  barrier  width  by 
the  formation  of  a  Ti-MoS2  mixture  layer  at  the  interface  is  thought  to  be  the 
dominant  mechanism  of  the  decreased  contact  resistance  after  laser  annealing  as 
Ti-doped  M0S2  becomes  more  conductive  due  to  increased  n-type  doping.  Fig.  6(d) 
shows  the  schematic  energy  band  diagram  for  the  M0S2  transistor  at  the  thermal 
equilibrium  and  after  laser  annealing  process.  If  we  assume  the  existence  of  a 
Schottky  barrier  between  metal  contacts  and  a  semiconductor,  thermionic  emission 
and/or  tunneling  will  allow  electron  transport  during  transistor  operation.  According 
to  Das  et  al.,  thermionic  emission  and  tunneling  though  Schottky  barrier  limit  the 
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charge  injection  during  off-  and  on-state  of  M0S2  transistors,  respectively.  The 
dominant  mechanism  that  reduced  contact  resistance  is  not  clear  as  the  gate  was 
floated  during  the  measurement  of  source-drain  current.  However,  the  significant 
increase  of  the  on-current  after  laser  annealing  in  Fig.  6(a)  suggests  that  laser 
annealing  increases  the  tunneling  current  across  the  Schottky  barrier.  Since  the 
tunneling  current  will  exponentially  increase  with  Nd°^,  where  Nd  is  carrier 
concentration,  the  decreased  contact  resistance  after  laser  annealing  is  probably  due  to 
the  reduced  tunneling  barrier  width,  which  could  be  resulted  from  the  increased 
carrier  concentration  at  the  interface.  To  confirm  the  results  and  further  investigation, 
we  will  examine  the  effect  of  the  laser  annealing  using  YFM  in  the  next  section. 
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Figure  6:  Device  performance  of  M0S2  TFTs  on  flexible  PEN  substrates,  (a)  Transfer 
characteristics  (Ids-Vgs)  for  the  flexible  M0S2  transistor  in  log-scale  (left)  and  linear-scale  (right)  before 
and  after  laser  annealing  at  V^s  =  1  V.  (b)  Output  characteristics  (Ids-V^s)  before  (dot)  and  after  (solid) 
laser  process,  (c)  Magnifications  (before  and  after)  for  the  low  Vds  range  from  0  to  0.1  V.  (d)  schematic 
energy-band  diagrams  of  a  multilayer  M0S2  transistors  with  a  Schottky  barrier  at  equilibrium  and  after 
laser  annealing  process. 


iv.  Y-function  method:  intrinsic  mobility 

The  field-effect  mobility  in  the  linear  regime  (peffjin)  and  the  po,  which  characterizes 
the  maximum  available  mobility  without  a  contact  factor  in  the  transistors,  were 
estimated  in  order  to  understand  the  laser  annealing  impact  on  the  contact  conditions 
or  M0S2  channel  layers  of  the  fabricated  devices.  Normally,  to  extract  the  contact 
resistance  of  devices,  the  transfer-line  method  (TLM)  is  employed.  However,  TLM 
requires  several  sets  of  transistors  with  various  channel  lengths  and  uniform  contacts. 
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Furthermore,  TLM  is  quite  demanding  and  difficult  to  apply  in  the  case  of  the  devices 
with  arbitrary  or  uncontrollable  shapes  and  very  small  size  channels.  Due  to  this 
limitation,  YFM  has  been  proposed  for  the  extraction  of  the  low-field  mobility 
without  the  influence  of  Rc  (i.e.  the  intrinsic  mobility,  po)  and  of  other  parameters  (e.g. 
Rc,  mobility  attenuator  factor  (0),  Vth)  as  well.  For  the  YFM,  there  is  no  need  to  make 
deliberated  and  calculated  devices;  the  YFM  can  prove  to  be  a  simple  and  powerful 
method.  The  YFM  is  based  on  the  Ids-Vgs  relation  in  the  linear  region,  as  below: 

W 

^ds  ~  8ds  ^  ^ds  ~  ~j^Qcht^eff  ^  ^ds  1  ' 

where  gds  is  the  channel  conductance,  peff  is  the  effective  mobility,  and  Qch  is  the 
channel  charge  per  unit  area.  Then,  we  consider  the  0  (=0o+0  ),  including  the 
contribution  from  the  channel  interface  (0o)  and  the  contact  resistance 
(0  =poCoxRcW/L).  Therefore,  the  Ids  can  be  written  as, 

W 

^ds  Sds  ^  ^ds  ~J~Qchl^eff  ^  ^ds 

Considering  the  definition  of  the  transconductance  ^^ds/dVgs,  Yds  =  const, 
the  Peff  (=  Lgm/(WCoxVds)),  the  Y  function  can  be  defined  as 


w 


~j~Coxtivyds  X  Y th) 


(5) 


in  order  to  eliminate  0.  Note  that  Eq.  (5)  is  independent  of  Rc,  which  is  assumed  to  be 
constant.  It  can  be  seen  from  Eq.  (5)  that  po  could  be  extracted  from  the  slope  of  the  Y 
function,  as  shown  in  Eig.  7(a).  Using  slopes  obtained  before  (5.9  x  10'"^)  and  after 
(6.6  X  10'"^)  laser  annealing,  the  values  of  po  are  extracted.  Eig.  7(b)  shows  the 
comparison  of  the  po  and  the  peak  peff  at  the  same  Vds,  +1.0  V.  Before  the  laser 
process,  a  large  discrepancy  (48.18%)  exists  between  po  (37.80  cm^  V'^s  ')  and  peffjin 
(19.59  cm^  V'^s'^).  However,  after  laser  annealing,  the  discrepancy  (2.95%)  is  greatly 
reduced:  po  (47.31  cm^  V's'^)  and  peffjin  (45.91  cm^  V'^s'^).  To  further  investigate  the 
laser  effects,  we  extracted  the  values  of  0,  and  estimated  Rc  by  the  following 
equations: 


Sffi  ~ 


dlo 


dYc 


Vp=cona 


-c 

L  [I  +  0{yc  —  Yt)^ 


X  Yd 


(6) 


Eor  obtaining  0,  Eq.  (6)  can  be  rewritten  as. 


y/Sm 


(7) 


where  Gm=  (W/L)poCox  is  the  transconductance  with  the  po.  Erom  Eq.  (7),  the  values 
of  0  for  as-fabricated  (1.0  x  10'^)  and  after  laser  annealing  (6.3  x  10'^)  can  be  obtained 
by  the  linear  fitted  slopes  (before:  ~  17  and  after:  ~  9.6)  of  the  l/gm^^-Vo  curve,  as 
shown  in  Eig.  7(c).  Through  the  known  values  of  0,  the  upper  bound  of  Rc  (Rc_max) 
can  be  calculated  when  we  assume  that  the  influence  of  0o  is  negligible.  Therefore,  the 
contact  effect  mainly  contributes  to  0,  and  extracted  values  of  Rc_max  are  28.9  kQ  and 
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14.6  kQ,  respectively,  before  and  after  the  laser  process.  Furthermore,  po  is  slightly 
increased  (20.10%)  before  and  after  the  laser  treatment.  Since  po  is  not  affected  by  the 
contact  factor,  one  of  the  possible  reasons  for  this  could  be  that  the  quality  of  the 
interface  is  enhanced.  The  traps  between  semiconductors  and  insulators  play  a  critical 
role  in  determining  the  device  performance,  and  the  reduced  interface  trap 
concentrations  were  confirmed  by  the  deceased  SS  values.  To  reconfirm  the  results, 
we  roughly  estimated  the  values  of  channel  resistance  (Rch)  by  using  the  relation 
Rch=L/(jUoWCox(VG-VTH))  before  and  after  laser  annealing.  The  obtained  values  of  Rch 
were  85.34  kQ  (before)  and  61.18  kQ  (after),  respectively.  The  results  of  decreasing 
Rch  by  the  laser  process  supports  and  agrees  well  with  the  enhanced  interface 
conditions  we  described  above.  Also,  through  the  comparison  with  Rc_max  and  Rch,  we 
have  evidence  that  the  diminishing  discrepancy  between  the  peffjin  and  po  originates 
from  the  reduced  contact  and  channel  resistance  by  the  laser  annealing  process.  This 
indicates  that  the  laser  annealing  process  employed  is  the  critical  factor  involved  in 
enhancing  the  TFT  performance  by  improving  the  contact  conditions  and  promoting 
the  interface  quality  without  thermal  damage  to  plastic  substrates. 


(b) 


Figure  7:  Y  function  and  evaluation  of  electrical  parameters  before  and  after  the  laser  process. 

(a)  Plots  of  Y  function  with  respect  to  Vgs  and  the  linear  fitting  in  the  strong  accumulation  region  for 
obtaining  the  slope,  (b)  Comparison  of  the  low-field  mobility  (po)  and  peak  field  effect  mobility 
(Peffjin)  before  and  after  the  laser  process,  (c)  Plots  of  with  respect  to  the  Vgs  and  their  slopes  in 

the  straight  line  fitting  for  extraction  of  the  mobility  attenuation  factor  (0). 


V.  Variability  of  electrical  contact  properties 

In  Korea  collaborators,  to  obtain  more  information  on  the  metal-MoSa  junction,  they 
measure  the  I-V  characteristics  of  the  two  M0S2  transistors  at  increased  temperatures 
up  to  240  °C.  Fig.  8(a)  shows  the  transfer  characteristics  of  the  two  M0S2  transistors 
at  five  different  temperatures  (room  temperature,  60  °C,  120  °C,  180  °C  and  240  °C). 
Since  the  two  transistors  have  different  channel  widths  and  thicknesses,  drain  current 
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density  (Jd)  is  used  for  comparison  instead  of  Id.  (Sample  A:  L  =  10  |J,m,  W  =  10  |J,m, 
and  thickness  of  M0S2  channel  t  =  31  nm,  Sample  B:  L  =  10  |im,  W  =  35  |J,m,  t  =  63 
nm). 

For  both  transistors,  off-current  density  (Joff)  and  subthreshold  swing 
significantly  increase  with  increased  temperature  as  Joff  is  known  to  be  dominated  by 
thermionic  emission  and  subthreshold  swing  is  proportional  to  kT  (k:  Boltzmann 
constant,  T:  temperature).  While  Joff  of  the  two  transistors  are  comparable  at  each 
temperature,  the  on-current  density  (Jon)  of  Sample  A  is  higher  than  that  of  Sample  B 
by  an  order  of  magnitude.  This  suggests  significant  difference  in  tunneling  barrier 
width  between  the  two  transistors  as  tunneling  is  known  to  dominate  the  on-current  of 
M0S2  transistors.  To  support  this  argument,  they  estimate  the  Schottky  barrier  height 
(Ob)  of  Samples  A  and  B  at  flat  band  condition.  They  build  a  series  of  Arrhenius  plot 
in  Figs.  8(b)  and  8(c)  with  the  data  in  Fig.  8(a).  Then,  assuming  the  conventional 
thermionic  emission  theory,  they  obtain  effective  Schottky  barrier  height  at  a  specific 
Vg  in  Fig.  8(d)  from  the  slope  of  the  data  in  Figs  8(b)  and  8(c).  In  thermionic  emission 
theory,  Jd  is  related  to  Ob  through  the  expression  Jd  =  A*  T^  exp(-qOB/kT) 
[exp(qVd/kT)  -  1],  where  A  is  Richardson  constant,  q  is  the  elementary  charge,  and  T 
is  temperature.  The  Ob  at  flat  band  condition  can  now  be  readily  estimated  in  Fig. 
8(d).  The  obtained  Ob  of  Samples  A  and  B  at  flat  band  condition  are  relatively  close 
(41  meV  and  50  meV,  respectively).  The  estimated  Ob  in  this  study  is  comparable 
with  the  reported  Ob  of  Ti-MoS2  in  literature  (50-65  meV). 


Figure  8:  (a)  Transfer  characteristics  of  two  M0S2  transistors  (Samples  A  and  B)  at  different 
temperatures.  Drain  current  density  (Jd)  is  used  for  comparison  due  to  the  differences  in  channel  width 
and  length  of  the  two  transistors,  (b)  Arrhenius-type  plot  of  Sample  A  created  by  the  data  in  (a),  (c) 
Arrhenius-type  plot  of  Sample  B  created  by  the  data  in  (a),  (d)  Extracted  effective  barrier  height  of 
Samples  A  and  B  as  a  function  of  applied  gate  voltage. 
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Fig.  9  shows  the  extracted  field  effect  mobility  of  Samples  A  and  B  from  room 
temperature  to  240°C.  At  this  temperature  range,  theoretical  calculation  indicates  that 
optical  phonon  scattering  dominates  the  transport  of  multilayer  M0S2  crystals  with 
ohmic  contacts.  The  mobility  of  Sample  A  slightly  decreases  from  31  cm^V's'^  to  29 
cm^V's'^  with  increasing  temperature  exhibiting  a  typical  signature  of  band  transport. 
However,  the  mobility  of  Sample  B  increases  exponentially  (from  5  cm^V'^s  '  to  11 
cm^V's'^)  for  the  same  temperature  range  suggesting  a  thermally  activated  transport. 
Three  possible  sources  can  explain  the  observed  thermally  activated  behavior  of 
Sample  B  in  Fig.  9:  the  existence  of  grain  boundaries  in  the  M0S2  channel,  the 
existence  of  Schottky  barriers  at  the  Ti-MoS2  junction,  or  the  existence  of  high 
concentration  of  trap  states  in  the  M0S2  channel.  With  the  existence  of  grain 
boundaries  ruled  out  due  to  the  single  crystal  M0S2  channel,  the  more  likely  source  is 
the  existence  of  Schottky  barriers.  Then,  the  distinct  transport  behavior  in  Fig.  9  can 
be  understood  by  the  energy  band  diagram  shown  in  the  inset  of  Fig.  9. 
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Figure  9:  Field  effect  mobility  of  Samples  A  and  B  as  a  function  of  inverse  temperature.  Inset  shows 
schematic  energy  band  diagrams  of  Samples  A  and  B  at  the  Ti-MoSi  junction. 


When  the  conduction  band  moves  upward  at  negative  gate  voltages,  the 
increased  barrier  height  prevents  electron  transport  from  metal  to  M0S2.  Since  this 
off-current  is  dominated  by  thermionic  emission  and  both  samples  have  similar 
barrier  height,  the  off-current  in  Fig.  9(a)  is  similar  for  each  sample.  When  the 
conduction  band  moves  down  at  positive  gate  voltages,  the  barrier  height  remains  the 
same  for  each  sample  due  to  Fermi  level  pinning.  Yet,  as  the  on-current  is  dominated 
by  thermally-assisted  tunneling.  Sample  A  will  have  a  much  narrower  barrier  width 
than  that  of  Sample  B  allowing  an  enhanced  tunneling  current  across  the  barrier. 
Therefore,  the  mobility  of  Sample  A  shows  the  dominant  contribution  of  M0S2 
channel  resistance  with  negligible  effect  of  the  Ti-MoS2  junction.  The  existence  of 
highly-doped  region  such  as  metallic  TiS2  or  Ti-doped  M0S2  has  been  observed  at  the 
Ti-MoS2  junction.  However,  due  to  the  wide  Schottky  barrier,  the  mobility  of  Sample 
B  is  dominated  by  the  contact  resistance  at  the  Ti-MoS2  junction. 

Since  the  tunneling  current  will  exponentially  increase  with  Nd°  where  Nd  is 
carrier  concentration,  significant  difference  in  carrier  concentration  will  exist  at  the 
Ti-MoS2  junction  between  Sample  A  and  Sample  B.  Such  a  variable  carrier 
concentration  at  the  Ti-MoS2  junction  could  be  related  to  the  non-uniform  distribution 
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of  metallic  defects,  assumed  to  be  clusters  of  S  vacancies,  in  naturally-occurring 
M0S2  crystals.  The  variation  of  defect  density  in  M0S2  could  induce  the  variation  of 
carrier  concentration  in  M0S2,  which  can  cause  the  variation  of  Schottky  barrier  width 
at  the  Ti-MoS2  junction.  Hence,  the  variation  of  defect  distribution  in  M0S2  crystals 
may  eventually  lead  to  the  different  behavior  of  mobility  depicted  in  Fig.  9. 

b.  Laser  doping 

i.  Laser  doping  of  monolayer  M0S2  and  WSe2:  tuning  the  optical  properties 
The  PL  spectrum  of  monolayer  M0S2  and  WSe2  before  and  after  the  laser  phosphorus 
doping  is  shown  in  Fig.  10(a).  The  peak  PL  intensity  of  the  laser-doped  area  is 
approximately  order  of  magnitude  greater  and  lesser  than  that  of  the  as-prepared 
monolayer  M0S2  and  WSe2,  respectively.  The  spectral  distribution  of  the  PL  peak  also 
sharpens  upon  the  laser  phosphorus  doping  while  the  peak  energy  of  the  PL  spectra  in 
monolayer  M0S2  was  blue-shifted  by  approximately  32  meV,  as  explained  by  the 
contributions  of  the  exciton  (X,  e-h  pair)  and  the  negative  trion  (X-,  e-e-h  clusters)  in 
the  PL  peak.  The  PL  peak  location  is  consistent  with  the  positions  of  the  exciton  peak 
(X;  -1.88  eV,  e-h  pair)  and  the  negative  trion  peak  (X-;  -1.84  eV,  e-e-h  pair).  For  the 
as-prepared  monolayer,  the  PL  contribution  of  the  negative  peak  (X-)  that  is 
positioned  at  1.84  eV  is  higher  than  that  of  the  exciton  peak  (X),  positioning. 
However,  in  the  PL  spectra  of  monolayer  M0S2  after  phosphorus  doping  by  laser,  the 
exciton  peak  (X)  contribution  dominates  the  effect  of  the  negative  trion  peak  (X-)  due 
to  recombination  of  excess  carriers  (electrons  in  M0S2)  by  hole  injection  from  the 
incorporated  phosphorus.  The  blue  shift  of  the  PL  spectrum  is  thus  consistent  with 
p-doping  effect.  The  p-doping  effect  on  intrinsically  n-type  M0S2  displays  largely 
different  effects  (Fig.  10(b))  compared  to  intrinsically  p-type  materials,  WSe2  (Fig. 
10(c)).  The  as-exfoliated  monolayer  WSe2  exhibits  direct  band  gaps  at  1.67  eV.  Upon 
phosphorus  doping,  the  PL  intensity  of  WSe2  is  drastically  reduced  since  incorporated 
excess  hole  carriers  strongly  suppress  the  exciton  emission.  The  spectral  shape  of 
peak  becomes  broader  after  the  laser  phosphorus  doping.  The  PL  spectra  peaks  in 
monolayer  WSe2  was  red- shifted  approximately  45  meV,  as  can  be  explained  by 
considering  the  contributions  of  the  exciton  (X,  e-h  pair)  and  the  positive  trion  (X-I-, 
e-h-h  clusters). 


Figure  10:  Photoluminescence  (PL)  spectra  of  monolayer  WSe2  and  M0S2  before  and  after  laser 
assisted  doping  with  phosphorus  as  a  p-type  dopant  (a)  and  its  PL  mapping  images  of  WSe2  (b)  and 
M0S2  (c). 
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The  Gaussian  profile  of  the  laser  beam  intensity  is  expected  to  generate  a 
spatially  varying  temperature  distribution  at  the  irradiation  spot.  Although  various 
reports  have  been  published  on  optical  and  thermal  properties  of  TMDC  materials, 
there  is  still  insufficient  information  to  enable  accurate  prediction  of  the  induced 
temperature  in  our  laser  doping  process.  The  maximum  PL  intensities  across  the  laser 
doped  regions  in  the  monolayers  M0S2  (Fig.  11(a)  inset)  and  WSe2  (Fig.  11(b)  inset) 
are  plotted  in  Figs.  11(a)  and  11(b).  These  distributions  could  be  fitted  to  a  Gaussian 
profile,  suggesting  that  the  hole  doping  concentration  increases  with  doping 
temperature,  which  maximizes  at  the  center  of  the  laser  reaction  spot.  The  full  width 
half  maximum  (FWHM)  of  the  PL  intensity  profile  was  approximately  2.5  pm,  which 
is  in  agreement  with  the  focal  spot  size  of  ~2.3  pm,  considering  the  IQx  and  0.28  NA 
lens  and  the  532  nm  laser  light  wavelength.  Fig.  11(c)  inset  shows  that  for  M0S2,  the 
PL  gradually  shifts  weight  from  X'  to  X  emission  as  the  reaction  spot  is  scanned  from 
the  edge  to  the  center.  In  contrast,  the  WSe2  shows  a  transition  from  X  to  X'^  (Fig. 
11(d)  inset). 


Spot  position 


Figure  11:  PL  peak  intensity  line  profile  of  laser  assisted  doping  monolayer  (a)  M0S2  and  (b)  WSea  at 
the  reaction  spot.  Inset  images  show  the  PL  peak  intensity-mapping  image  and  optical  microscopy 
image  of  the  laser-doped  region.  Scale  bars  are  2  and  5  pm,  respectively.  Variation  of  the  PL  intensity 
within  the  single  spot  of  laser  doped  monolayer  (c)  M0S2  and  (d)  WSe2.  Inset  shows  the  PL  mapping 
image.  The  spacing  between  each  spot  is  500  nm. 


a.  Nano-Auger  electron  spectroscopy 

To  study  the  mechanisms  of  the  laser-assisted  doping  process,  nano-Auger  electron 
spectroscopy  (nano-AES)  which  enables  local  chemical  analysis  using  focused 
electron  beam  of  ~10  nm  spot  size  at  ultra-high  vacuum  of  base  pressure  -10  '°  mbar 
(Oxford/Omicron  system)  was  performed  on  pristine  and  laser-doped  regions  of  M0S2 
crystals  at  various  doping  conditions  such  as  irradiation  time,  doping  gas  composition. 
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and  laser  power.  Considering  that  AES  reflects  in  the  chemical  properties  of  the 
material  surface  within  a  depth  of  few  atomic  layers,  multilayer  and  tri-layer  M0S2 
flakes  were  chosen  as  model  strictures  analogous  to  monolayers  to  explore  the  role  of 
laser  heating  during  doping  process.  In  Fig.  12(a),  the  S/Mo  Auger  signal  ratio  of  the 
multilayer  M0S2  layer  is  slightly  reduced  by  10  sec  and  30  sec  irradiation  implying 
that  chalcogen  vacancies  are  created  at  low  density,  while  it  can  be  clearly  shown  that 
the  S/Mo  ratio  decreases  significantly  and  the  Auger  transitions  of  P  atoms  are  shown 
up  after  60  sec  irradiation.  The  drastic  changes  in  Auger  spectra  occurred  between  30 
sec  and  60  sec.  Therefore,  it  is  presumable  that  the  phosphorous  atoms  were  directly 
deposited  on  the  M0S2  surface  after  30  sec  laser  irradiation  via  the  dissociation  of  PH3 
dopant  molecules  accelerated  by  the  heat  accumulation.  Surface  roughness  change 
after  irradiation  for  60  sec  was  also  confirmed  by  AFM  scanning.  Also,  Fig.  12(b) 
clearly  shows  that  as  the  concentration  of  PH3  dopant  gas  balanced  with  Ar  is  higher 
(i.e.,  Ar  flow  rate  is  set  lower),  the  S/Mo  ratio  of  the  multilayer  M0S2  flake  is 
decreased  more  after  laser  irradiation. 

Further,  the  laser-assisted  doping  reaction  in  the  2D  regime  has  been  analyzed 
quantitatively  by  performing  nano-AES  on  M0S2  tri-layers  as  summarized  in  Fig. 
12(c).  The  P  doping  level,  [P]/([P]+[S]),  was  estimated  to  be  4.68  at.%  and  56.04 
at.%  from  the  tri-layers  that  were  laser-  doped  at  180  mW  and  200  mW  respectively 
while  the  doping  duration  time  was  set  as  60  sec  for  the  both  cases.  In  the  case  of  200 
mW  process,  P  atoms  seem  to  be  deposited  on  the  M0S2  surface  excessively,  which 
gives  a  very  similar  result  with  the  cases  when  the  laser  power  or  gas  composition  are 
set  higher  than  threshold  values  as  described  in  Figs.  12(a)  and  12(b). 
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Figure  12:  (a)  Auger  electron  spectroscopy  spectra  taken  before  and  after  the  doping  at  a  range  of  laser 
irradiation  times  (10,  30,  and  60  sec  at  180  mW  with  multilayer  M0S2),  (b)  dopant  concentration  by 
argon  dilution  (flow  rate  of  100  and  500  seem  with  multilayer  M0S2),  and  (c)  laser  irradiation  power 
(180  and  200  mW  at  60  sec  with  tri-layer  M0S2). 


Hi.  Laser  doping  of  bilayer  M0S2  and  WSe2:  tuning  the  electrical  properties 
In  parallel  to  strong  modification  of  the  luminescent  properties,  the  laser-assisted 
doping  also  imposes  drastic  effects  on  electrical  properties  of  the  ultra-thin  TMDCs. 
This  was  investigated  using  field  effect  transistor  (FET)  devices  incorporating  the 
locally-doped  TMDCs  layers  as  current  channels.  As  displayed  in  Figs.  13(a)  and 
13(c),  multiple  electrodes  was  patterned  onto  TMDC  bilayer  flakes,  covering  both 
pristine  and  laser-doped  regions.  A  voltage  applied  through  the  back  gate  provides 
additional  modulation  of  carrier  density  in  the  TMDC  channel.  In  the  undoped  regions, 
n-type  (p-type)  FET  behavior  was  observed  with  on-off  ratio  exceeding  five  orders  of 
magnitude  from  the  M0S2  (WSe2)  bilayers.  The  FET  characteristics,  i.e.,  source-drain 
current  (Isd)  versus  back  gate  voltage  (Vg)  curves,  of  devices  fabricated  with  different 
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distance  from  the  center  of  the  laser-doped  region,  change  gradually  in  a  manner 
consistent  with  the  spatially-graded  variation  of  PL  across  the  same  region,  as 
summarized  in  Figs.  13(b)  and  13(d).  The  phosphorous  doping  suppresses  the  n-type 
characteristics  of  the  M0S2,  and  enhances  the  p-type  characteristics  of  the  WSea 
drastically,  with  increase  in  conductance  as  well  as  shift  of  threshold  voltage.  The 
hysteresis  features  of  FET  loops  are  likely  due  to  water  molecule  absorption  on  the 
surface  of  the  TMDCs,  considering  that  all  electrical  measurements  were  conducted 
under  ambient  conditions. 
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Figure  13:  Bilayer  M0S2  and  WSe2FET  device  (a),  (c)  the  flakes  are  outlined  by  the  blue  dashed  box 
and  the  laser  irradiated  area  is  denoted  by  the  red  dashed  circle.  Source-drain  current  (Ld)  versus  back 
gate  voltage  (Vg)  curve  recorded  from  the  device  in  (a)  and  (c)  with  a  bias  voltage  of  SOOmV  (b),  (d) 


iv.  Long-term  stability  of  the  laser-assisted  phosphorus  doped  M0S2 
Reliable  performance  of  devices  requires  stable  doping  of  the  active  materials.  Indeed, 
previous  studies  have  already  achieved  both  n-  and  p-type  doping  by  surface  charge 
transfer  mechanisms,  mostly  through  chemical  physisorption.  However,  physisorption 
doping  is  unstable  and  decays  almost  completely  within  an  hour,  or  is  retained  for 
longer  periods  only  if  the  doping  environment  is  maintained  or  protected.  Among  the 
various  doping  methods,  substitutional  doping  of  foreign  elements  is  an  effective  and 
stable  doping  strategy  for  TMDCs.  In  our  work,  the  laser  was  used  to  create  sulphur 
vacancies  and  locally  heat  up  the  material  to  crack  the  precursor  molecules,  enabling 
the  substitution  of  sulphur  with  phosphorous.  The  PL  data  in  Fig.  14(a)  shows  that  the 
laser-assisted  doping  is  irreversible  and  stable,  even  after  exposure  in  air  for  a  month. 
As  a  comparison  to  test  the  doping  stability,  we  also  studied  a  laser  annealing  effect. 
In  this  case,  the  monolayer  M0S2  was  laser  irradiated  in  absence  of  the  dopant  gas  at 
the  given  laser  conditions  (i.e.  the  exact  same  conditions  used  for  doping).  As  shown 
in  Fig.  14(b),  the  PL  intensity  of  monolayer  M0S2  increases  right  after  the  laser 
annealing,  as  in  the  case  of  laser  doping.  However,  the  PL  intensity  increase  was 
higher  than  the  increase  in  the  case  of  laser  doping.  After  one  week,  however,  the  PL 
intensity  drastically  drops  and  continues  to  drop  to  the  original  value  within  2  weeks, 
indicating  instability  of  the  laser  annealing  effect,  and  the  necessity  of  the  presence  of 
the  PH3  molecules  for  stable  doping. 
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Figure  14:  (a)  PL  spectra  of  laser  assisted  phosphorus-doped  monolayer  M0S2  and  (b)  laser  annealed 
with  the  same  condition  of  doped  one.  The  samples  are  under  exposure  to  oxygen  and  water  molecules 
in  ambient  air. 


c.  Laser  welding 

i.  Plasmonic  laser  welding 

The  electrical  contact  at  the  Ag  NW  junctions  can  be  improved  by  laser  welding, 
which  is  based  on  highly  localized  heating  induced  by  plasmonic  effects.  Fig.  15(a) 
illustrates  the  schematic  concept  of  the  laser  welding  process  for  the  2D-like  random 
networks  of  Ag  NWs.  In  this  research,  as-prepared  2D-like  Ag  NW  networks  showed 
unfavorable  electrical  sheet  resistance,  480  ~  550  O/n,  due  to  weak  physical  contact 
between  nanowires  (in  the  left-bottom  inset  of  Fig.  15(a)).  Here,  we  note  that  even 
though  different  sheet  resistance  can  be  produced  by  tuning  the  nanowire  diameter, 
length,  and  surface  coverage  of  nanowires,  we  applied  laser  process  under  ambient 
condition  to  enhance  the  electrical  properties  of  the  interconnected  wires  (in  the 
right-bottom  inset  of  Fig.  15(a)),  without  requiring  precise  control  of  the  dimensions 
or  the  number  density  of  the  Ag  NWs. 

After  the  laser  welding  process,  the  sheet  resistance  was  drastically  reduced  by 
10  ~  15  Q/n;  these  values  were  almost  37  or  55-  fold  lower  than  for  the  as-prepared 
samples.  To  confirm  this  finding,  we  took  SEM  images  before  and  after  laser 
irradiation  as  shown  in  Figs.  15(b)  and  15(c).  The  results  clearly  indicate  that  the  laser 
process  induced  fusion  (red  circles)  at  the  NW  junctions  (Fig.  15(c))  in  contrast  to  the 
as-prepared  NW  junctions  (Fig.  15(b)).  The  features  of  laser  welding  process  for 
creating  the  interconnected  mesh  are  well  described  by  considering  the  plasmonic 
effect  at  Ag  NW  junctions  that  form  a  naturally  small  gap.  The  concentrated  photon 
source  induces  localized  heat  spreads  into  the  NWs  due  to  their  high  thermal 
conductivity.  The  small  gaps  at  the  NW  junctions  enable  selective  and  effective 
focusing  of  light  and  heating  (to  sufficiently  high  temperature  to  soften/melt  the  outer 
NW  surfaces),  right  at  the  contact.  We  could  find  evidence  of  the  locally  generated 
plasmonic  effect  generated  in  the  small  gap  configuration  from  the  morphological 
changes  of  the  NWs.  As  shown  in  the  inset  of  Fig.  15(c),  after  the  laser  illumination, 
the  morphology  of  only  the  NW  junction  substantially  changed;  NWs  slightly  melted 
and  joined  via  welding.  On  the  other  hand,  the  morphology  of  NWs  away  from  the 
junction  did  not  change  unlike  in  the  junction  proximity.  If  the  thermal  distribution 
were  uniform  on  the  NWs,  morphological  changes  would  be  observed  over  their 
entire  lengths.  Therefore,  this  laser  welding  process  was  different  from  conventional, 
uniformly  distributed  thermal  heating.  The  locally  generated  plasmonic  effect  can  also 
prevent  damage  to  low  thermal  budget  plastic  substrates. 
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Figure  15:  Plasmonic  nanowelded  metal  NWs  network  for  flexible  gate  electrode,  (a)  The 

schematic  concept  of  the  laser  welding  process  for  Ag  NW  network  on  the  PI  substrate  formed  by 
filtration  method.  The  bottom  insets  show  schematic  illustrations  of  the  junctions  between  two  adjacent 
NWs  [Left]  before  and  [Right]  after  laser  treatment,  (b-c)  SEM  images  of  (a)  as-prepared  Ag  NWs 
random  networks  before  laser  welding  and  (c)  after  the  laser  treatment  (red  arrows  indicate  the  laser 
welded  junctions).  The  inset  shows  the  magnified  SEM  image  of  laser- welded  part. 


a.  Solution  based  embedded  PI  structure 

Laser  welding  treatment  makes  conductivities  of  the  Ag  NW  films  less  susceptible  to 
external  mechanical  stress,  but  considerable  surface  roughness  still  persists.  A 
possible  solution  to  this  issue  is  to  embed  the  Ag  NW  electrode  into  the  substrate, 
which  also  eliminates  the  film-substrate  adhesion  problem. 

For  a  flexible  substrate,  solution-type  PI  is  spin-coated  onto  the  laser-treated 
Ag  NWs  network,  resulting  in  the  permeation  of  PI  into  the  voids  formed  between  the 
Ag  NWs  owing  to  its  low  viscosity  and  surface  energy.  As  a  result,  the  Ag  NWs  are 
partially  embedded  and  held  firmly  within  the  flexible  PI  substrate  with  their  initial 
positions  and  corresponding  contacts  remaining  intact  under  varying  mechanical 
stress. 

Fig.  16(a)  compares  the  atomic  force  microscope  (AFM)  images  of 
as-deposited,  normal  (i.e.,  not-embedded)  Ag  NWs  on  a  rigid  glass  substrate  and 
partially  embedded  Ag  NWs  in  a  PI  substrate.  Individual  NWs  are  clearly  identified  in 
the  top  AFM  image  in  Fig.  16(a).  However,  the  surface  morphology  and  roughness 
change  dramatically  when  the  Ag  NWs  are  partially  embedded  in  the  PI,  exhibiting  a 
much  smoother  surface  with  reduced  roughness  as  shown  in  the  bottom  AFM  image. 
These  results  are  also  confirmed  in  the  topological  line  profile  (Fig.  16(b)). 
Root-mean-squared  heights  for  normal  and  embedded  Ag  NWs  are  found  to  be  ~39.0 
nm  and  6.0  nm,  respectively. 


DISTRIBUTION  A.  Approved  for  public  release:  distribution  unlimited. 


Figure  16:  Solution  based  embedded  polyimide  (PI)  structure,  (a)  2D  AFM  images  of  the  Ag  NWs 
as-deposited  either  on  a  rigid  glass  substrate,  i.e.,  not  embedded  (top),  or  partially  embedded  in  the  PI 
substrate  (bottom),  (b)  Topological  line  profile  data  obtained  along  the  black/blue  dashed  lines  in  (a). 


Hi.  Numerical  simulation  of  nonlinear  hyperelastic  materials 

The  mechanical  behavior  of  hyperelastic  materials  (e.g.,  PI,  PEN,  PET,  polymer,  or 
rubber-like  materials)  is  generally  expressed  in  terms  of  strain  energy  potential. 
Eurthermore,  many  different  models  have  been  suggested  to  define  the  strain  energy 
function  for  predicting  the  properties  of  hyperelastic  materials.  Therefore,  in  order  to 
come  up  with  a  proper  model  that  describes  well  the  mechanical  characteristics  of  the 
flexible  M0S2  devices  with  hyperelastic  substrate,  different  hyperelastic  models  must 
be  compared  with  the  real  experimental  data.  Eig.  17(a)  shows  the  measured 
stress(cr)-stretch(A)  curve  of  a  specimen  with  a  thickness  (t)  of  0.105  mm,  fabricated 
with  solution  based  PI  and  PET.  Under  static  and  cyclic  bending  test,  the  strain  (a) 
represents  the  amount  of  force  per  unit  area,  and  the  engineering  stretch  (A)  represents 
the  relative  elongation: 


^  (8) 

where  s  is  the  engineering  normal  strain.  Erom  the  experimental  data,  an  ultimate 
tensile  stress  of  46.6  x  10  N  m-  is  observed  at  a  maximum  elongation  of  43%  (A  = 
1.43).  Also,  within  the  moderate  strain  range  of  interest  (<  10%,  A  =  1.10),  the 
Mooney-Rivlin  model  is  found  to  match  well  with  our  experimental  data  as  shown  in 
the  inset  of  Eig.  17(a).  The  two-parameter  Mooney-Rivlin  model  assumes  that  the 
local  strain  energy  density  in  an  incompressible  material  is  a  simple  function  of  local 
strain  invariants.  Eor  uniaxial  loading,  the  nominal  stress-stretch  behavior  for  the 
Mooney-Rivin  model,  is  given  by: 

where  Aq  is  the  cross-section  area  of  the  test  specimen  and  Lq  is  the  original  length  of 
the  specimen.  The  two-parameters,  Cio  and  Coi,  are  material  parameters  determined 
by  fitting  the  model  (Eq.  (9))  to  the  experimental  data  (blue  dot  in  the  inset  of  Eig. 
17(a)).  Through  the  model,  the  stress  state  in  a  cross  section  of  the  flexible  substrate 


DISTRIBUTION  A.  Approved  for  public  release:  distribution  unlimited. 


mounted  on  a  mechanical  bending  tester  can  be  simulated,  as  shown  in  Fig.  17(b). 
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Figure  17:  Computational  non-linear  hyperelastic  analyses  on  the  mechanical  bending  stress  of 
the  flexible  M0S2  TFT,  (a)  Measured  stress((T)-stretch(l)  curve  of  the  specimen  of  solution  based  PI 
and  PET  (t  =  0.105  mm).  The  inset  shows  experimental  data  fitted  with  Mooney-Rivlin  model,  (b)  The 
numerical  result  on  hyperelastic  nonlinear  stress  analyses  based  on  the  parameters  extracted  from  the 
experimental  data  when  the  bending  radius  is  5  mm.  The  top  inset  shows  a  photograph  of  the  bent 
M0S2  TFTs  with  a  radius  of  5  mm.  (c-d)  Based  on  the  Mooney-Rivlin  model,  the  expected 
mechanically  stable  range  (green  region)  of  the  structure  with  respect  to  the  (c)  radius  of  curvature  of 
the  substrate  and  (d)  thickness  of  the  substrate. 


To  investigate  the  flexible  structure,  hyperelastic  nonlinear  stress  analysis  is 
performed  through  COMSOL  Multiphysics.  It  is  generally  known  that  when  an 
elastomeric  slab  is  bent  with  arbitrary  radius  r,  there  exists  a  neutral  plane  {i.e., 
stress-free  zone)  in  the  geometric  mid-plane.  The  computational  result  shows  that  the 
outer  surface  of  the  specimen  is  under  tensile  stress  and  the  inner  surface  in  under 
compression,  with  respect  to  the  neutral  plane  where  the  material  experiences  no 
stress  (Fig.  17(b)).  When  the  bending  radius  (r)  was  5  mm,  the  maximum  internal 
tensile  stress  induced  in  the  outermost  region  (interface  between  Ag  NWs  network 
and  PI  substrate)  is  4.38  x  10^  N  m-^.  This  induces  a  total  strain  of  0.2%  in  the 
substrate.  Such  small  strain  induced  in  the  PI  upon  bending  indicates  that  the  devices 
also  undergo  relatively  small  mechanical  strain  levels,  and  so  the  electrical  properties 
are  affected  to  a  very  small  extent.  To  understand  the  correlation  between  the  bending 
radius  and  the  induced  stress  in  detail,  calculations  were  performed  with  respect  to  the 
bending  radius  and  the  substrate  thickness.  Through  such  a  procedure,  one  may 
predict  a  range  of  device  dimensions  over  which  the  system  is  mechanically  stable.  As 
show  in  Fig.  17(c),  the  flexible  structure  is  anticipated  to  tolerate  a  specific  range  of 
stress  (green  region,  r  >  2.5  mm  at  t  =  0.105  mm)  without  reaching  the  failure  stress 
(4.66  X  10^  N  m-^)  and  elongation  (43%).  Note  that  the  results  based  on  the 
Mooney-Rivlin  model  do  not  guarantee  device  stability  at  relatively  large  strains 
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(marked  with  dotted  lines).  Also,  when  the  substrate  thickness  t  is  larger  than  0.8  mm 
(at  r  =  5  mm),  the  structure  can  be  destroyed  (Fig.  17(d)).  Therefore,  in  order  to 
remain  in  the  mechanically  stable  state,  the  generated  stress  must  be  kept  below  the 
ultimate  tensile  strength  (red  region).  Additional  parameters,  such  as  the  thicknesses 
and  mechanically  properties  of  each  individual  layer  {e.g.,  dielectric,  semiconductor, 
electrodes,  etc.),  as  well  as  their  structure  must  also  be  optimized  so  that  functional 
devices  can  be  achieved.  When  ductile  materials  such  as  metals  having  very  low 
elastic  region  (only  0.2  or  0.3%  strain)  are  used,  the  structure  becomes  more 
vulnerable  to  the  bending  so  a  more  meticulous  design  is  mandatory.  In  this  regard, 
laser  welded  Ag  NWs  network  can  thus  be  preferential  over  metal  layers,  not  only  for 
the  flexible  but  also  stretchable  electronics. 


iv.  Variation  of  TFT  performances  under  mechanical  stresses 

For  investigating  solid  tolerance  of  our  flexible  M0S2  TFTs  against  various 
mechanical  stresses,  two  kinds  of  tensile  stresses  were  applied  to  the  flexible  TFTs 
array  parallel  to  their  active  channels;  static  and  cyclic  upward  bending  (Insets  of  Figs. 
18(a)  and  18(b)).  The  lamination  of  a  100-pm-thickness  PET  film  beneath  the  PI 
substrate  enables  the  downshift  of  neutral  plane  in  the  flexible  M0S2  TFTs  array,  so 
that  the  entire  TFT  could  lie  within  the  tensile-stress  zone.  As  a  result,  the  flexibility 
of  our  devices  can  be  investigated  under  relatively  harsher  environments. 


Figure  18:  Variation  of  M0S2  TFT  performances  under  various  mechanical  stresses,  (a) 

Comparison  of  the  transfer  characteristics  of  the  flexible  M0S2  TFTs  under  static  bending 
environments  with  three  different  bending  radii  (flat,  r  =  10  mm,  and  5  mm).  The  inset  shows  a 
photograph  of  the  bent  M0S2  TFTs  attached  on  rigid  cylinder  with  a  pre-defined  radius  (white  dashed 
arc  and  arrow),  (b)  Comparison  of  the  transfer  characteristics  under  cyclic  bending  with  respect  to 
number  of  bending  (n  =  0,  20,  100,  and  1,000).  The  inset  shows  a  photograph  of  the  flexible  M0S2 
TFTs  loaded  on  a  multi-modal  bending  tester.  The  cyclic  bending  radius  was  fixed  at  10  mm  (black 
dashed  circle  and  arrow),  (c-d)  Variations  of  /Jefi  and  Vh  of  the  flexible  M0S2  TFTs  as  a  function  of  the 


DISTRIBUTION  A.  Approved  for  public  release:  distribution  unlimited. 


(c)  static  bending  radii  and  (d)  number  of  cyclic  bending. 

Fig.  18(a)  compares  the  transfer  characteristic  curves  of  the  M0S2  TFTs  array 
under  flat  and  statically  bent  conditions.  The  hs-Vgs  curves  of  the  M0S2  TFT  with 
bending  radii  of  10  and  5  mm  did  not  change  significantly  with  respect  to  that  in  the 
flat  condition.  The  variation  of  the  mobility  values  (A/t  (%)  =  l/t  bending  -  /tflatl  //fflat)  due 
to  the  bending  with  r  =  10  and  5  mm  are  estimated  to  be  2.0%  and  8.8%,  respectively 
(Fig.  18(c)).  The  threshold  voltage  (Vth)  of  the  M0S2  TFT  in  the  flat  condition  is  -18.7 
V.  For  r  =  10  mm,  the  Vth  just  remains  identical  to  its  initial  value,  and  shifts  to  -19.5 
V  when  r  =  5  mm.  Fig.  18(b)  shows  a  comparison  of  the  transfer  characteristics  of 
another  M0S2  TFT  under  cyclic  bending.  As  the  number  of  bending  (n)  increases,  the 
entire  hs-Vgs  curves  shift  slightly  toward  negative  Vgs  value.  However,  the  hn  /  7off 
ratio  is  preserved  to  ~10^.  The  Aju  (Vh)  due  to  the  cyclic  bending  with  n  of  20,  100, 
and  1,000  are  1.9%  (-10.1  V),  5.3%  (-11.3  V),  and  10.0%  (-15.4  V),  respectively  (Fig. 
18(d)).  It  should  be  noted  that  the  thin  ALD  AI2O3  layer  plays  an  important  role  in  our 
flexible  TFTs  array,  by  enhancing  the  adhesion  between  the  GI  and  the  S/D  electrodes 
as  well  as  active  layer.  M0S2  TFTs  without  AI2O3  layer  exhibited  poor  transistor 
performance  under  identical  bending  conditions. 


d.  High-mobility  transistors  based  on  large-area  and  highly  crystalline  CVD-grown  2D  films 

The  uniform  growth  of  large  grain  multilayer  MoSe2  films  on  insulator  substrates  was 
demonstrated  using  a  modified  chemical  vapor  deposition  method.  Highly  crystalline  MoSe2 
film  with  large  grain  sizes  was  fabricated  as  an  active  channel  for  a  transistor,  consisting  of 
bottom  gate  architecture  with  Ti/Au  source  and  drain,  on  flexible  plastic  substrate  or  a  rigid 


Figure  19:  Synthesis  and 


microstructure  of  MoSca  film  and  electrical  characteristics  of  a  flexible  MoSea 
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transistor  under  mechanical  bending  test. 

Our  multilayer  MoSe2  field-effect  transistors  exhibited  high  mobility  (~  121  cm  V 
's~  ^),  a  high  on/off  current  ratio  (>  lO"^),  and  good  mechanical  flexibility  while  bending  down 
to  radii  of  5  mm.  The  experimental  results  can  be  understood  remarkably  well  by  classical 
long-channel  MOSFET  theory.  The  MoSe2  films  presented  in  this  work  are  attractive  for  the 
fabrication  of  high-speed  thin-film  devices  for  flexible  integrated  circuits,  especially  for 
applications  such  as  humanoid  robotics,  flexible  OLED  displays,  and  human-centric  soft 
electronics.  This  paper  was  published  in  journal  of  ‘Advanced  Materials’,  which  was  selected 
as  a  cover  image. 

e.  Real-time  electrical  detection  of  epidermal  skin  M0S2  biosensor  for  point-of-care 
diagnostics 

Toward  point-of-care  diagnostics,  various  approaches  have  been  proposed  and,  in  particular 
optical  detection  is  preferred  due  to  its  relatively  simple  and  rapid  process.  At  the  same  time, 
EET-based  biosensors  have  attracted  great  attention  as  they  can  provide  highly  sensitive  and 
label-free  detection.  In  this  work,  we  present  highly- sensitive  point-of-care  devices  in  an 
epidermal  skin-type  with  system-level  integration  of  flexible  M0S2  EET  biosensor,  read-out 
circuits,  and  light-emitting  diode  (LED)  that  allow  real-time  detection  of  prostate  cancer 
antigen  (PSA).  Regardless  of  physical  forms  or  mechanical  stress  conditions,  our  proposed 
high-performance  M0S2  biosensors  can  detect  PSA  concentration  of  1  pg/ml  without  specific 
surface  treatment  for  anti-PSA  immobilization  on  the  M0S2  surface.  Eurthermore,  current 
modulation  induced  by  the  binding  process  was  stably  maintained  longer  than  2~3  minutes. 
The  demonstrated  results  indicate  that  flexible  MoS2-based  EET  biosensors  have  a  great 
potential  as  point-of-care  diagnostics  of  prostate  cancer  as  well  as  other  biomarkers.  This 
manuscript  was  selected  as  a  cover  image  of  Nano  Research  (2017,  volume  4). 


Figure  20:  The  cover  image  and  an  illustration  and  demonstration  of  platform  of  an  epidermal  skin-type  M0S2 
biosensor  system:  (a)  Schematic  of  a  flexible  M0S2  biosensor  detecting  label-free  immunoassay  of  PSA,  where 
PSA  antigen  directly  binds  with  PSA  antibody  functionalized  M0S2  surface,  (b)  Schematic  layout  of  2D 
multilayer  M0S2  FETs.  c,  d)  Schematic  illustration  (c)  and  photograph  (d)  of  an  epidermal  skin-type  M0S2 
biosensor  system,  consisting  of  biosensor,  read-out  circuits,  and  LED  as  an  indicator. 


DISTRIBUTION  A.  Approved  for  public  release:  distribution  unlimited. 


In  a  first  demonstration,  M0S2  biosensors  are  designed  to  electrically  detect  PSA  in  a 
highly  sensitive  and  label-free  manner  without  the  need  of  a  chemically  treated  gate 
dielectric.  The  nature  of  hydrophobic  M0S2  surface  (the  contact  angle  ~  75.77)  allows  a 
physical  adsorption  to  biomolecule  binding,  and  doesn’t  require  additional  dielectric  layers 
such  as  Hf02  or  Si02.  Such  an  oxide-free  operation  can  provide  higher  device  sensitivity  as 
well  as  a  simpler  device  structure.  Moreover,  the  use  of  “OFF”  current  for  a  transistor  as  an 
indicator  strongly  affords  additional  improvement  in  device  sensitivity  in  M0S2  biosensor  and 
can  compute  the  surface  charge  concentration  due  to  anti-PSA/PSA  binding  quantitatively 
(Figs.  21(a)  and  21(b)).  Theoretical  calculation  of  the  M0S2  biosensor  operation  provides  a 
systematic  understanding  on  the  experimental  demonstration  of  the  M0S2  biosensors  in  this 
study,  and  will  open  new  pathways  to  fundamental  studies  of  electrical  2D  sensor 
applications  for  biological  and  chemical  species. 
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Figure  21:  (a)  Schematic  of  a  M0S2  biosensor  configured  as  a  PSA,  detecting  label-free  immunoassay, 
illustrating  PSA  antibody  functionalized  M0S2  surface  (top)  and  subsequent  binding  of  PSA  antigen  with 
antibody  receptors,  (b)  Change  of  the  off-current  versus  various  PSA  concentrations  for  an  anti-PSA  modified 
n-type  M0S2  transistor.  As  the  concentration  of  PSA  increases  from  1  pg/mL  to  10  ng/mL,  the  amount  of 
selective  binding  of  PSA  to  the  anti-PSA  immobilized  M0S2  nanosheet  is  approximately  proportional  to  the 
PSA  concentration. 
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